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Abstract 

The perturbative QCD corrections to the semileptonic decay width of the 
tau lepton are evaluated in the next-next-to-leading order from the contour 
integral representation in various renormalization schemes, using numerical 
solution of the renormalization group equation for complex energies. A quan- 
titative estimate of the ambiguities resulting from the freedom of choice of 
the renormalization scheme is obtained by taking into account predictions in 
all schemes that do not involve large cancellations in the expression for the 
scheme invariant combination of the expansion coefficients. The problem of an 
optimal choice of the renormalization scheme for the improved perturbative 
expression is discussed. A fit of is made using the available experimental 
data. 
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As is well known, predictions of a quantum field theory obtained with finite 
order perturbative expressions depend on the choice of the renormalization scheme 
(RS). This effect is usually regarded as a complication in the comparison of the 
theory with the experimental data. However, the effect of renormalization scheme 
dependence may be also used to improve our understanding of the perturbative 
result. By estimating the change in the predictions over a set of a priori acceptable 
schemes we may test reliability of the perturbative expression. This point of view 
has been recently discussed in H, where a concrete method has been proposed 
for a systematic evaluation of the RS ambiguities in the QCD predictions in the 
next-next-to- leading order (NNLO). This method is based on the observation, that 
one should take into account predictions in all schemes that do not involve large 
cancellations in the expression for scheme independent combination of the expansion 
coefficients. Application of this method was illustrated in ^ using as an example 
the QCD corrections to the Bjorken sum rule for the polarized structure functions. 



Also the problem of the RS ambiguities in the QCD corrections to the ratio []3|-||10||: 
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has been considered in this framework [|I|]. In this case it was found, that even small 
changes in the RS - among acceptable schemes - resulted in a large variation in the 
predictions. The strong RS dependence seemed to undermine the phenomenological 
importance of the QCD corrections to Rr- 

However, the QCD correction to R^ comes originally in the form of a contour 
integral in the complex energy plane Pj-f^l- The commonly used perturbative for- 
mula appears as a result of approximate evaluation of the contour integral. As was 
pointed out in [O, O, the evaluation of the contour integral may be improved. In 



it was shown that using the improved evaluation of the contour integral one ob- 
tains smaller scale dependence in the MS scheme. It is then an interesting problem, 
how the improved accuracy of the predictions affects the estimate of the strength of 
the RS dependence evaluated according to the method developed in [||, |]. 

In this note we evaluate the QCD corrections to Rr in various renormalization 
schemes by solving numerically the renormalization group equation for complex 
energies, and by computing numerically the relevant contour integral. Using the 
method developed in fl], |2| we obtain a quantitative estimate of the strength of the 



RS dependence in this case [[T^. We discuss the problem of an optimal choice of 
the renormalization scheme in the improved expression. We also make a fit of A^g^ 
taking into account recently obtained experimental results for i?^. 

Let us begin with a brief summary of the ideas presented in H, 0]. Let us consider 
a NNLO expression for a physical quantity 5, properly normalized, depending on a 
single energy variable P: 

5{P^) = a(p2)[i + r,a{P^) + r^a'iP'')], (2) 
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where a(yU^) = (y'^(yU^)/(47r^) denotes the running couphng constant that satisfies the 
NNLO renormahzation group equation: 

/i— = — (1 + CiO + C2a^ ), (3) 
a/i 

(The effects of the quark masses are neglected.) The value of the predictions ob- 
tained for 6 with the NNLO expressions (@) and (H) depends on the choice of RS. 
Presently in the phenomenological applications the modified minimal subtraction 
scheme (MS) |]I4[ is used. However, other choices of the renormahzation scheme are 
possible, all related to the MS scheme by a finite renormahzation. In the NNLO, in 
the class of the mass and gauge independent schemes, the freedom of choice of the 
scheme may be characterized by two continuous parameters. (The relevant formulas 
describing the RS dependence of the coefficients ri, r2 and C2 have been collected 
in [^.) The consistency of the perturbation theory guarantees, of course, that dif- 
ferences in the predictions obtained in various schemes are formally of higher order 
in the coupling constant, but numerical magnitude of these differences is significant 
for phenomenological analysis. This fact stimulated the interest in prescriptions for 



making a proper choice of the RS [15|-|24]. The most attractive of these propositions 



seems to be the choice based on the so called Principle of Minimal Sensitivity (PMS) 
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However, whatever scheme we decide to choose as an optimal one, there is a con- 
tinuum of equally reasonable schemes close to the one preferred by us. Predictions 
in such schemes also should be somehow taken into account in the phenomenological 
analysis. A natural way to do this is to supplement the prediction in a preferred 
scheme with an estimate of the variation of the predictions over the whole set of 
a priori acceptable schemes. A concrete realization of this idea was presented in 
[0], based on the existence of an RS independent combination of the expansion 
coefficients [[T^, gT|, ^, |2|: 



p2 = C2 + r2 - Ciri - r^. (4) 

which provides a natural RS independent characterization of the magnitude of the 
NNLO corrections for the considered physical quantity. In it was proposed to 
calculate variation of the predictions for 6 over the set of schemes for which the 
expansion coefficients satisfy the condition: 

0"2(ri,r2,C2) < / IP2I, (5) 

where 

(J2{ri,r2,C2) = \c2\ + \r2\ +Ci|ri| + r^. (6) 

A motivation for the condition (^ is that it eliminates schemes in which the ex- 
pressions (0) and (^) involve unnaturally large expansion coefficients, implying large 
cancellations in the expression for the RS invariant p2- The constant / in the condi- 
tion (IP controls the degree of cancellations that we want to allow in the expression 
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for p2- As was pointed out in [Q], taking / = 2 in one has the PMS scheme right 
at the boundary of the allowed region in the space of parameters characterizing the 
scheme. 

It is expected that the estimate of the strength of the RS dependence obtained 
according to this prescription would be useful for a quantitative comparison of re- 
liabihty of perturbative predictions for different physical quantities, evaluated at 
different energies. It should be also very useful in determining the regions of appli- 
cability of the perturbative expansion. Indeed, any large variation of the predictions 
over a set of schemes satisfying constraint (j^) with / = 2 would be an unambiguous 
sign of a limited applicability of the NNLO expression. 

Let us now consider the QCD corrections to the semileptonic decay of the tau 
lepton. These corrections, denoted further by 5^, contribute to the ratio in the 
following way: 



Rr = mudV + \VusV)Sew{1 + ^EW + ^r) 



(7) 



1.0194, 
26]. The 



where Vij are the Cabbibo-Kobayashi-Maskawa matrix elements and Sew 
Sew = 0.001 are factors arising from the electroweak corrections | 
perturbative correction to 6r is known up to NNLO 1^-0 in the approximation of 
massless u, d and s quarks, and it has in the conventional approach the form (H) 
with the characteristic energy scale P = = 1.7771 GeV. 



The interest in the QCD corrections to Rr 0-||T2[, ||2^-[^ was originally stim- 
ulated by the fact, that the perturbative expression for 6r appeared to be quite 
sensitive to the QCD parameter A. Thus in principle a very good test of perturba- 
tive QCD predictions could be obtained ^ |^, since the effects of nonzero masses 
of light quarks, and the nonperturbative effects, were estimated to be small despite 
the rather low energy scale. Unfortunately, in [|l| large differences in the perturba- 
tive predictions for Sr have between obtained in a priori acceptable schemes. (The 
perturbative QCD corrections to Rr in some optimal schemes were considered pre- 
viously in [^, |28|, ^ |3^-) It seemed therefore that the high sensitivity of 6r to the 
QCD parameter A is practically eliminated by the ambiguities introduced by the 
strong RS dependence 



However, the commonly used perturbation expansion for Sr involves not only the 
obvious approximation resulting from the truncation of the perturbation series, but 
also an approximation of another kind, which appears in the process of evaluation of 
S-r from the expression for the two-point quark current correlators. Indeed, in order 
to evaluate the QCD corrections to Rr one starts from the formula ||, |], 0: 
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where n*^^-* denotes the sum of the transverse parts of the AS* = 0,1 vector and 
axial quark current correlators: 



n«(s) = |Kdr[n, 



(1) 

udy 



s)] + \vM^%[ 



^) + nl^i»], 



(9) 
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(q) = i-g^-q' + q'-qn^^ly/Aq' 



Klv/Aiq) = ^ / d'^xe^'^ < o\TiJj:,y^^ix)j;:,y/Aoy)\o > . 



(10) 

(11) 



(If the quark mass effects are neglected, the longitudinal part of 11'^'' does not con- 
tribute. Also, the electroweak contributions have been neglected.) Since the QCD 
predictions for the quark current correlators are not very well known for real posi- 
tive s, one uses the analyticity properties of 11 to convert the expression (||) into the 
contour integral in the complex energy plane: 
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where C is a contour running clockwise from s = rn^ — ie to s = m^ + ie and avoiding 
the region of small i. e. on C we have |s| = m^. Integrating by parts along the 
contour, one obtains: 



6r = — r de (l + 2e'^ 
111 J -It ^ 
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where (5n(— s) is defined via the relation: 
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-127r^)s-n5;)(s) = 3(|Kd|' + |K/)[1 + 5n{-s) 



(13) 



(14) 



The importance of this expression lies in the fact, that 5r has been expressed via the 
quantity (5n(— s), which is formally RS independent, and which is directly calculable 
in perturbative QCD for real negative s. (5n(— s) has a perturbation expansion of 
the form (|^), with the expansion coefficients in the MS scheme, with /i^ = — s, given 
by @: _ _ 

rf ^ = 1.63982, rf ^ = 6.37101. (15) 

The commonly used expansion of 5r in terms of a{'m'^) is obtained from the expression 
(13) by expanding a{—s) under the integral in terms of a{m1): 
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(16) 



and performing the contour integral over the appearing powers of ln(— s/m^) explic- 
itly. In this way one obtains the NNLO expression for 6r which has the form (H) 
with the expansion coefficients ff. 



ri 
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(17) 
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In the MS scheme we have 

= 5.2023, ff^ = 26.366. (19) 
We see that there is a significant difference between the expansion coefficients and 

It should be noted however, that using the expansion of a(— s) in terms of a(m^) 
in the NNLO expression for 6n, we make a rather crude approximation. For example, 
we effectively assume, that: 

a\-s) ^ a^ml). (20) 

As was pointed out in |jTl], |l^ one may obtain a more accurate description of the 
QCD effects in the tau lepton decay by evaluating the s-dependence of the running 
coupling constant on the contour C in a more precise way. 

It is then an interesting question, whether the improved evaluation of the QCD 
corrections may change the rather upsetting result found in [0. 

In this note we analyze the predictions for 6r obtained from the contour integral 
representation (|l^) using a numerical solution of the renormalization group equation 
in the complex energy plane. To obtain the running coupling constant along the 
contour C we solve numerically the implicit equation: 

= rp-ri + $(a,C2), (21) 




MS 



where 



$(a, ca) = ci In ( ) + - + Ci In(cia) + 0(a). (22) 
\2ciJ a 



The explicit form of $(a,C2) is given for example in |43]. (The renormalization 
group coefficients for nj = 3 are b = 4.5, ci = 16/9 and C2^^ = 3863/864 fa 4.471.) 
This equation is obtained by integrating the renormalization group equation (^) 
with an appropriate boundary condition |T^, and then analytically continuing to 



general complex s. (The presence of rf^"^ and in this general expression results 
from taking explicitly into account the one-loop relation between A parameters in 
different schemes |jl5|, which is valid to all orders of perturbation expansion. All 
values of Ap^ mentioned in this paper correspond to three "active" quark flavors.) 
To obtain the prediction for 6r we take under the contour integral the quantity Sn 
in the NNLO approximation in the form (0) and we perform the contour integral 
numerically. To analyze the RS dependence of the thus obtained expression for Sr 
we parametrize the freedom of choice of the RS by the coefficient ri in and the 
coefficient C2, and we calculate variation in the predictions when these parameters 
are changed in a region determined by the condition (|]). The region of allowed 
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values for ri and C2 for p2 > which is the case needed here, was described 

analytically in ||^. 

In our analysis a fundamental role is played by the NNLO RS invariant, which 
for 6yi has the value p2 = 5.238. This value should be compared with p2 = —5.476, 
which is obtained for the conventional expansion of 6r in terms of a(m^). The fact 
that p2 has similar magnitude, but the opposite sign, compared to p2 indicates that 
the approximations used to derive the conventional expansion for 6r may greatly 
distort the pattern of the RS dependence of the predictions. 

The results of our calculations are presented in four figures. In Fig.l we compare 
perturbative QCD predictions for Sr as a function of m^/A;^, obtained in the next- 
to-leading order (NLO) and in NNLO using the conventional expansion, with the 
predictions obtained by evaluating numerically the contour integral. We clearly see 
that the differences between NLO and NNLO are much smaller for the predictions 
obtained from numerical calculation of the contour integral than for the predictions 
obtained with the conventional expansion in terms of a(m^). This confirms the 
expectation that the procedure for evaluating the 6r adopted in this note gives 
indeed improvement over the conventional expansion. 

In Fig. 2 we show a contour plot for 6r as a function of ri and C2, obtained for 
^JTs ~ 310 MeV, together with the region of scheme parameters satisfying the 
condition with 1 = 2 — smaller region — and I = 3. (Note that — for purely 
technical reasons — the actual variable on the vertical axis is C2 — CiVi instead of 
C2.) We see that the pattern of RS dependence in this case has more complicated 
structure than that obtained for conventional NNLO approximants — this is evi- 
dent if one compares Fig. 2 with the corresponding figure in [^, representing the RS 
dependence of the QCD corrections to the Bjorken sum rule. In particular, we find 
here three critical points which are relatively close to the allowed region, instead of 
one encountered in the conventional approximant. The coordinates of the critical 
point closest to the origin, i.e. the saddle point at the boundary of the / = 2 allowed 
region, are very close to the PMS parameters for the quantity 6u evaluated for real 
negative s ("euclidean" PMS). (Let us note that because of a nonpolynomial char- 
acter of the improved expression for 6r the value of the PMS parameters cannot be 



obtained from the set of algebraic equations given in |17], but it has to be deter- 



mined from direct numerical analysis.) This confirms the fundamental importance 
of Su for 6r. It also shows that approximations used to derive the conventional 
NNLO expression greatly distort the pattern of the RS dependence — the conven- 
tional NNLO approximant for 5^ would have a critical point for completely different 
values of the scheme parameters. To estimate the strength of the RS dependence we 
compare predictions for 6r obtained for the scheme parameters lying in the denoted 
allowed region. Because of the location of the critical points in this case the maxi- 
mal and minimal values for Sr are attained at the boundaries of the allowed regions 
corresponding to / = 2 and / = 3. It is interesting to note that the MS scheme with 
C2 — ciTi ^1.56 lies slightly outside the / = 3 allowed region. 

In Fig. 3 we show, how the pattern of the RS dependence of the improved ex- 
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pression for 5r depends on the value of Xnimr / ly-jjg) ■ We see that for small values 
of we have only one critical point in the / = 2 allowed region, and that it lies 
very close to the PMS parameters for 5n evaluated for spacelike momenta. With 
increasing Ajg^ the structure of the critical points close to the / = 2 allowed region 
becomes more complicated. For some values of three critical points appear. 
For large values of ^-j^ the critical point corresponding to the "euclidean" PMS 
evolves into a very flat plateau, and the secondary critical point, with negative ri, 
moves into the / = 2 allowed region. 

In Fig. 4, which is the main result of this paper, we show how the optimal per- 
turbative predictions for 5^-, and the estimates of the RS dependence, behave as a 
function of rrir /A.jjg. The results of our calculations are compared with constraints 
from the experimental data. As our preferred predition, indicated in Fig.4 by a thick 
solid line, we take the values obtained in the scheme with ri = 0, C2 = 1.5p2- These 
parameters describe, for small and intermediate values of A;^, the approximate 
location of the leading critical point in the (^1,02) plane, and for larger values of 
A;^ they fall in the region of very small variation of the predictions. From Fig. 3 it 
is clear, that for all interesting values of Ajjg the predictions for 6r in this scheme 
practically coincide with the values at the relevant critical points. The dashed and 
dash-dotted lines represent the variation of the predictions over the I = 2 and I = 3 
allowed regions, respectively. In order to make our results useful for other authors 
we give in Table 1 the numerical values of the preferred predictions for 6r, together 
with an estimate of the RS dependence, for several values of \n{mT- / Ajj^) . For com- 
pleteness we include also the NLO predictions, evaluated using the contour integral 
and optimized according the the PMS prescription — they are well approximated 
by taking the scheme with ri = —0.76. 

We see that the variation of the predictions over the / = 2 allowed region has 
a rather smooth dependence on A;^, being almost constant for the considered 
range of Ajjg. The variation over the / = 3 region of parameters is reasonably 
close to the I = 2 variation, except for the maximal value, which grows rapidly for 
larger A;^. This should be compared with the large differences found in [|l| for the 
conventional NNLO approximant for 5^. We see that using an improved evaluation 
of the predictions for 6r we qualitatively improve stability of the predictions with 
respect to change of the RS, and we greatly reduce the ambiguities resulting from 
the freedom of choice of the renormalization scheme. 



Let us now compare our results with those obtained in ||T2[. The analysis of 
[0 concentrates on variation of the renormalization scale /i in the MS scheme, 
which in our approach corresponds to variation of ri for fixed C2 = c^^ ■ In 
only a brief comment is made on the variation of predictions when C2 is changed 
from to 2c2^^ with fixed fi. It must be stressed, that there is no theoretical or 
phenomenological reason justifying a priori such a choice of the set of the scheme 
parameters. Also, there is no immediate relation between the magnitude of variation 
found according to [|1^] and the estimate of the RS dependence obtained according 
to our method. In particular, small RS dependence obtained by prescription of 
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Q does not necessarily imply a small variation in our approach. It should be 



emphasized, that if we agree that there is a "democracy" of renormalization schemes, 



then the set of considered schemes must be larger than that discussed in , and any 
reasonable constraint on the freedom of choice must be related to the RS invariant p2- 
Most importantly, a consideration of an essentially two-dimensional set of scheme 
parameters is unavoidable if one intends to find the optimal NNLO predictions 
according to the Principle of Minimal Sensitivity, which requires for the optimal 
scheme the vanishing of the derivatives with respect to both scheme parameters. 

The results presented in Fig. 4 may be used to obtain an improved constraint 
on Ajj^ from the world average of the available experimental data on Rr- As is 
well known, the experimental value of Rr may be obtained in three independent 
ways: from the electronic and muonic branching ratios in the tau lepton decay, 
and from the total width of the tau lepton. (Note that Particle Data Group gives 
two sets of values for the leptonic branching fractions — the number shown below 
corresponds to the set called "our average." See Appendix for a detailed discussion 
of the subtleties involved in obtaining experimental value for Rr.) Taking a weighted 
average of these three determinations we find: 

{RrTeJp = 3.591 ± 0.036. (23) 

This value may then be converted into experimental constraint on 6r via Eq.(^. 
However, in order to have a phenomeno logically meaningful constraint on (5^^'"*^ we 
must take into account also the effect of the light quark masses and the nonper- 



turbative effects. According to |T^, ^ these may be estimated to give roughly an 
overall —(1.5 ± 0.4)% correction to the value of Rr- Therefore, after taking into 
account that the factor involving the CKM matrix elements is practically equal to 
unity, our final formula relating Rr and 5^^''* takes the form: 

Rr = 3x 1.0194 X (1.001 + 5f *■*) X (0.985 ± 0.004). (24) 

Using ( P^D we obtain: 

(5r'^*)exp = 0.191 ± 0.012(exp) ± 0.005(npt), (25) 

where the first error reflects the effect of experimental uncertainties and the second 
one the effect of uncertainties in nonperturbative corrections. 
Fitting Ai^ to this experimental value we obtain: 

Ag^ = 376(opt)l j^(th,l=2) ± 29(exp) ± 12(npt) MeV. (26) 

For the / = 3 region of parameters the variation in A^^ is MeV. 

As was explained in ^, a natural way to parametrize the QCD predictions 
in our analysis is via the A;^. However, for comparison with results given in 
other papers on experimental tests of QCD, we give also the corresponding values 
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of as(m^) in the MS scheme, obtained with the three-loop renormahzation group 
equation: 

a^{ml) = 0.332l°:°°^(th,l=2) ± 0.015(exp) ± 0.006(npt). (27) 

For the / = 3 region of parameters we obtain for af^^^ml) the variation of Iq oio- 

Using the formulas given in [^] to match Ap^ for different number of flavors we 
obtain the corresponding result for as(m|): 

= 0.1190t[!:[J[]E]^(th,l=2) ± 0.0017(exp) ± 0.0007(npt). (28) 

For the / = 3 region of parameters we obtain for a*^'^(m|) the variation of taoois- 
We see that the RS dependence ambiguities in the determination of a^^^^m"^) 
or a'f^{m\) from are not very big — this is the result of stabilization of the 
predictions when improved evaluation procedure is used. Note however that the 
RS dependence ambiguities are still comparable in magnitude to the uncertainties 
related to the present experimental accuracy of Rr- 

Let us note that using for the leptonic branching fractions the values denoted by 
PDG as "our fit" we obtain (see Appendix): 



{RrTeJp = 3.535 ± 0.037. (29) 
which is 1.5(7 smaller than (p3]). Using this value we find in NNLO: 



A^ = 331MeV(opt), af ^(m^) = 0.309, af ^(m|) = 0.1162, (30) 

with the same theoretical and experimental uncertainties as before. 

To have a complete picture of theoretical uncertainties we give also the result 
of a fit of A;^ to the experimental value of Rr, obtained with the improved NLO 
predictions optimized according to the PMS prescription. Using the NLO expression 
with ri = —0.76 we obtain: 

A^ = 413(opt) ± 32(exp) MeV. (31) 

Using then the two-loop renormahzation group equation we find: 

af^iml) = 0.316 ± 0.013(exp), (32) 

af^(m|) = 0.1209 ± 0.0018(exp). (33) 

Summarizing, we have evaluated the perturbative QCD corrections to the semilep- 
tonic decay width of the tau lepton, using the contour integral representation and 
evaluating numerically the running coupling constant on the contour. We obtained 
results in a broad class of renormahzation schemes. We discussed the problem of 
obtaining the optimal predictions in the improved evaluation, and we found predic- 
tions for St in the scheme preferred by the Principle of Minimal Sensitivity, both 
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in NLO and in NNLO. Using a specific condition to eliminate the schemes with 
unnaturally large expansion coefficients we obtained a quantitative estimate of the 
ambiguities in 5t- arising from the freedom of choice of the RS. We expect that this 
estimate would be useful in combining the QCD results from the Rr with other 
experimental constraints on QCD. We found that the improved expression for 5^- 
is much more stable with respect to change of RS compared to the conventional 
perturbative NNLO expression. Nevertheless, we found that the RS dependence 
ambiguities in the fit of A^g^ to the experimental data for 5r are at least of the order 
of uncertainties arising from the present accuracy of experimental determinations of 



Note added 

Soon after this paper was posted on the bulletin board, two papers appeared 



4^ in which determination of from the tau decay is discussed in detail. We make 



two comments about these papers. First, the value of obtained in is almost 
identical to the value given in our paper. However, this value was obtained using 
(i?T-)exp = 3.56 ±0.03, which is la lower than the value we have used in our analysis. 
Secondly, the main source of the theoretical uncertaintly estimated in ap- 
pears to be the contribution from the yet uncalculated 0(Q;f) corrections. Such an 
estimate necessarily involves some speculative assumptions on the magnitude of the 
higher order terms. Also, the higher order contributions would be RS dependent, 
with the freedom of choice of the RS characterized by three arbitrary parameters in 
the four-loop order. It is straightforward to extend the method used in our paper 
to the four-loop case in order to obtain a quantitative estimate of the resulting RS 
ambiguity. 
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Appendix 

In this appendix we make some remarks on the experimental determination of Rr- 
First, let us note that the experimental value of Rr may be obtained from three 
independent measurements: from the electronic and muonic branching ratios in the 
tau lepton decay, and from the total width of the tau lepton. To obtain R^ from 
the electronic branching ratio we take: 

{^)exp = 777^; ^ ~ \ 1^ ) ' 

\-De)exp \ J->e / th 
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where according to formulas given in with shghtly updated parameters, we 
should take: 

[-^\ = 0.972568. 

\BeJ th 

To obtain from the muonic branching ratio we take: 

exp (^B^)exp \BeJth V -^e 

Finally, to obtain Rr from the total decay width we take: 

/J-, _ (Xtot)exp _ ^ _ (^Ji 

with 



th 



{Te)th = (0.40341 ± 0.00057) x 10^'^GeV 



obtained from the formulas given in ||25[ with updated parameters. 

As may be expected, these three determinations do not give exactly the same 
value. Therefore we take for R =< R > ±6R a weighted average, according to the 
prescription: 

where Wi = {6Ri)~^. 

Secondly, one has to be careful which set of the world averaged experimental data 
one uses in the phenomenological analysis. The Particle Data Group [44] gives in fact 
two sets of values for the leptonic branching fractions. One set is a straightforward 
world average of the data: 

{Be)exp = 0.1790 ± 0.0017, {B^)exp = 0.1744 ± 0.0023, 

This set is called by PDG "our average." Another set, 

iBe)exp = 0.1801 ± 0.0018, (5^)exp = 0.1765 ± 0.0024, 

called "our fit," is a result of a global fit to the tau lepton branching fractions, taking 
into account the relevant constraints. Taking "our average" of PDG for the leptonic 
branching fractions we obtain: 

(^r)fxp = 3.614 ± 0.053, {Rr)X = 3.604 ± 0.074. 

Taking the present world average for the total decay width: 

{Ttot)exp = = (2.227 ± 0.023) x lO-^^Q^y, 

yTrJexp 



13 



we find: 

(it:,)^;^^ 3.548 ±0.065. 

(For Tr the PDG also gives two values — in this case it seems however clear that 
one should take the value in which results of older experiments have been corrected 
for the change in the experimental number for m^.) Taking a weighted average of 
three values for we find: 

{Rr)Zl = 3.591 ± 0.036. 

This value has been further used in our paper to obtain constraints on A^g^. 

It has to be emphasized however, that if we use "our fit" of PDG for the leptonic 
branching fractions, we find: 

{Rr)exp = 3.525 ± 0.055, {Rr)% = 3.538 ± 0.075. 

Together with the value obtained from the total decay width this gives a weighted 
average of: 

{Rr)Zl = 3.535 ± 0.037. 
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Figure Captions 



Fig.l. Perturbative QCD predictions for 5r as a function of rrir /Ajj^, obtained 
in NLO and NNLO with the conventional expansion (dashed hnes) and with the 
numerical evaluation of the contour integral (solid lines). 

Fig. 2. Perturbative prediction for 6r, obtained from the numerical evaluation of the 

(3) 

contour integral for =310 MeV, as a function of the parameters ri and C2 
determining the RS. Note that the actual variable on the vertical axis is C2 — ciVi. 
The regions of scheme parameters satisfying the condition (^ with / = 2 and / = 3 
are also indicated. At the critical point inside the / = 2 region we have 6r = 0.1649. 



Fig. 3. Same as in Fig. 2, but for various values of ln(m^/A^g^). 

Fig. 4. Perturbative prediction for 6r, obtained from the numerical evaluation of the 
contour integral representation, as a function of m^/A^^. The thick lines denote the 
preferred NNLO predictions (i.e. obtained in the scheme with ri = and C2 = 1.5p2 
— see text) and the optimal NLO predictions (obtained with ri = —0.76). Also 
indicated are the maximal and minimal values obtained in NNLO when the scheme 
parameters are varied over the / = 2 (dashed lines) and / = 3 (dash-dotted hnes) 
allowed regions. For comparison we show the experimental constraint on the value 
of (5f 



15 





{^t)nNLO 


Xrnin 
^1=2 


Xmax 
"l=2 


Xmin 
^1=3 


Xmax 
"l=3 


('^-r)iVLO 


1.30 


.23634 


.2293 


.2443 


.2075 




.21816 


1.35 


.22652 


.2194 


.2322 


.2090 




.20909 


1.40 


.21711 


.2101 


.2223 


.2052 




.20056 


1.45 


.20816 


.2014 


.2135 


.1969 


.2401 


.19256 


1.50 


.19968 


.1932 


.2050 


.1890 


.2209 


.18505 


1.55 


.19168 


.1856 


.1971 


.1817 


.2055 


.17800 


1.60 


.18414 


.1785 


.1895 


.1749 


.1928 


.17139 


1.65 


.17705 


.1718 


.1823 


.1685 


.1853 


.16518 


1.70 


.17039 


.1656 


.1754 


.1625 


.1783 


.15935 


1.75 


.16413 


.1596 


.1689 


.1568 


.1718 


.15387 


1.80 


.15824 


.1541 


.1628 


.1515 


.1655 


.14871 


1.85 


.15271 


.1489 


.1569 


.1465 


.1596 


.14386 


1.90 


.14751 


.1440 


.1514 


.1418 


.1540 


.13928 



Table 1: Numerical values of the preferred NNLO predictions for 6r (i.e. obtained 
in the scheme with ri — and C2 = 1.5p2 — see text), and the maximal and 
minimal values obtained after variation of the scheme parameters in the allowed 

(3) 

region corresponding to / = 2 and / = 3, for several values of ln(m^/A^g^). For 
comparison also the values of the preferred NLO predictions are given (ri = —0.76). 
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